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Preface 

The  work  described  herein  was  performed  during  the  period  June  1985 
through  January  1986  by  Dr.  Ronald  A.  Erchul  under  an  Interagency  Personnel 
Agreement  (IPA)  with  the  Virginia  Military  Institute,  Lexington,  Virginia. 

The  work  was  funded  by  the  Repair,  Evaluation,  Maintenance,  and  Rehabilita¬ 
tion  CREMR)  Research  program  sponsored  by  Headquarters,  US  Army  Corps  of  Engi¬ 
neers  (USAGE) .  The  work  was  performed  under  Civil  Works  Research  Work 
Unit  32315,  "Geophysical  Techniques  for  Assessment  of  Existing  Structural 
Foundations."  The  REMR  Overview  Committee  consists  of  Mr.  James  E.  Crews  and 
Dr.  Tony  C.  Liu,  Headquarters,  USAGE.  The  Technlal  Monitor  was  Mr.  Ben  Kelly. 

Field  work  was  performed  by  Dr.  Erchul  and  students  of  his  engineering 
geology  classes  at  the  Virginia  Military  Institute.  This  report  was  prepared 
by  Dr.  Erchul.  Mr.  William  F.  McCleese  is  REMR  Program  Manager,  and 
Mr.  Jerry  S.  Huie  is  the  problem  area  leader.  Principal  Investigator  for  this 
REMR  research  project  is  Dr.  Dwaln  K.  Butler,  Earthquake  Engineering  and  Geo¬ 
physics  Division  (EEGD) ,  Geotechnical  Laboratory  (GL) ,  US  Army  Engineer  Water¬ 
ways  Experiment  Station  (WES).  Dr.  Arley  G.  Franklin  is  Chief,  EEGD,  and 
Dr.  William  F.  Marcuson  III  is  Chief,  GL.  The  report  was  edited  by 
Ms.  Odell  F,  Allen,  Information  Products  Division,  Information  Technology 
Laboratory,  WES. 

COL  Dwayne  G.  Lee,  CE,  is  Commander  and  Director  of  WES.  Dr.  Robert  W. 
Whalln  is  Technical  Director. 
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Conversion  Factors,  Non-SI  to  SI  (Metric) 


Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  (met¬ 
ric)  units  as  follows: 


Multipl 


To  Obtain 


degrees  (angle) 

0.017A5329 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins* 

K 

feet 

0.3048 

metres 

> 

KV 

inches 

2.54 

centimetres 

1 

miles  (US  Statute) 

1.609347 

kilometres 

*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  the  following  formula:  C  =  (5/9)  (F  -  32).  To  obtain  Kelvin  (K)  read 
ings,  use:  K  =  (5/9)  (F  -  32)  +  273.15. 
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GEOTECHNICAL  APPLICATIONS  OF  THE 
SELF  POTENTIAL  (SP)  METHOD 


THE  USE  OF  SELF  POTENTIAI.  IN  THE  DETECTION  OF  SUBSURFACE 
FLOW  PATTERNS  IN  AND  AROUND  SINKHOLES 


Introduction 


1.  The  purpose  of  this  research  was  to  test  and  evaluate  an  Indirect 
geophysical  measuring  technique  called  self  or  spontaneous  potential  (SP)  on  a 
newly  developed  sinkhole  to  determine  If  this  method  would  be  able  to  detect 
the  direction  of  ground-water  flow.  This  study  represents  a  novel  application 
of  the  SP  method;  and  this  report,  preliminary  In  nature,  presents  a  proof 
test  of  the  application  and  concept. 

2.  The  SP  method  measures  natural  electrical  potential  field  differ¬ 
ences  at  the  earth's  surface.  Anomalies  In  the  electric  field  can  be  produced 
by  ore  bodies,  heat  flow,  or  subsurface  fluid  flow.  SP  has  been  used  widely 
In  oil  well  logging  to  detect  permeable  zones,  and  It  has  also  been  used  to 
delineate  conductive  ore  deposit.  The  SP  method  has  been  used  for  many  years 
In  the  USSR  for  geotechnical  applications  such  as  seepage  analysis  and  study 
of  landslide  processes.  Publication  of  papers  by  Soviet  researchers  Ogllvy, 
Ayed,  and  Bogoslovsky  (1969)  and  Bogoslovsky  and  Ogllvy  (1972)  led  to  a  number 
of  geotechnical  applications  In  the  United  States. 

3.  Recently  the  US  Army  Engineer  Waterways  Experiment  Station  (WES)  has 
used  SP  measuring  techniques  at  Gathright  Dam  In  Virginia  (Cooper  and  Koester 
1982),  Clearwater  Dam  in  Missouri  (Koester  et  al.  1984)  and  Mill  Creek  Dam  In 
Walla  Walla,  Washington,  (Butler,  Wahl,  and  Sharp  1984)  to  determine  the 
abnormalities  created  In  the  ambient  electric  field  In  cases  where  water  was 
known  to  be  flowing  through  porous  zones  of  subsurface  material  but  the  loca¬ 
tion  of  these  zones  was  not  well  defined.  WES  has  greatly  simplified  the  SP 
measuring  technique  by  using  copper  clad  steel  grounding  rods  as  electrodes 
rather  than  the  nonpolarizing  electrodes  that  many  Investigators  have  used  in 
the  past.  For  the  scope  of  the  Army  projects  the  copper  clad  electrodes  were 
very  cost  effective,  and  they  detected  relative  potential  difference 
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successfully.  It  was  this  simplification  of  the  SP  measuring  technique  that 
encouraged  the  present  investigation. 

4.  The  application  of  SP  in  this  study  is  based  on  the  established 
observation  that  the  amount  and  the  direction  of  fluid  (ground-water)  flow 
sets  up  electromotive  forces  in  the  surrounding  media  with  resulting  varia¬ 
tions  in  potential  which  are  detectable  as  anomalies  at  the  surface.  Negative 
anomalies  are  associated  with  seepage  areas  (Ogilvy,  Ayed,  and  Bogoslovsky 
1969)  and  with  water  flowing  in  solution  channels  (Cooper  and  Koester  1982). 
Water  flowing  in  a  solution  channel  at  a  depth  of  90  ft  in  Manatee  Springs, 
Florida,  was  readily  detectable  by  SP  (Cooper  1983).  In  Cooper  and  Koester 
(1982)  negative  anomalies  of  500  mv  were  observed  and  assumed  to  be  paths  of 
ground-water  flow.  Electrical  resistivity,  dye  injections,  and  bore  holes 
confirmed  this  model. 

5.  When  water  flows  through  underground  channels,  the  flowing  fluid 
forms  a  double  boundary  layer  at  the  channel  wall  which  results  in  a  net  nega¬ 
tive  electrical  charge  in  the  geological  material  surrounding  the  channel 
(Butler,  Wahl,  and  Sharp  1984).  This  accumulation  of  negative  ions  at  the 
solid-liquid  phase  interface  is  then  detected  as  anomalous  to  the  positively 
charged  surrounding  rock  bodies,  and  these  areas  are  taken  to  be  the  preferred 
paths  followed  by  ground  water  as  it  flows  through  conduits,  fissures,  and 
fractures. 

Site  Description 

6.  The  study  area  is  6  miles*  south  of  Lexington,  Virginia,  in  the 
southern  Shenandoah  Valley.  This  is  an  area  of  karst  topography  and  contains 
many  sinkholes,  which  range  in  diameter  from  8  to  300  ft  and  in  depth  from 

3  to  50  ft.  Within  the  study  area  of  2  miles  by  1/2  mile,  25  sinkholes  occur 
(see  Figure  1).  Two  of  these  sinkholes  were  instrumented  for  this  study 
(Harris  and  Hunter  sinkholes) .  The  Harris  sinkhole  had  experienced  a  cata¬ 
strophic  collapse  during  a  heavy  rainfall  in  August  1984;  hence,  it  displays  a 
fresh  aven  (opening)  at  its  bottom.  The  Hunter  sinkhole  is  a  shallow  feature 
typical  of  subsidence  development. 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 
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Figure  1.  Topographic  map  of  the  study  area  and  the 
sinkholes  In  this  area 


7.  Developed  in  the  upper  portion  of  the  Ordovician  aged  Beekmantown 
Formation,  the  Hunter  and  Harris  sinkholes  occur  at  an  elevation  of  approxi¬ 
mately  1,200  ft.  The  Beekmantown,  a  thick  bedded,  finely  grained  dolomite  is 
approximately  1,500  ft  chick  in  the  study  area.  The  structure  of  the  area  is 
in  the  form  of  an  anticlinal  fold  that  plunges  to  the  southwest  (Spencer 
1968) .  The  Beekmantown  is  considered  to  be  one  of  the  most  Important  cave 
bearing  strata  in  Virginia.  Some  of  the  more  notable  commercial  caves,  (End¬ 
less  and  Skyline)  as  well  as  many  of  Virginia’s  long  caves,  occur  in  the  Beek 
mantown;  within  a  mile  of  the  study  area  there  are  five  caves. 

8.  The  area  is  mostly  open  pasture,  but  exposed  bedrock  outcrops 
(karren)  of  the  Beekmantown  Formation  are  common  in  the  study  area.  These 
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exhibit  evidence  of  extensive  solutional  weathering,  especially  along  joints. 
Base  level  is  Buffalo  Creek,  a  perennial  stream  located  1/2  mile  northeast  of 
the  study  area  at  an  elevation  of  approximately  1,000  ft  (200  ft  below  the 
studied  sinkholes) . 

9.  The  two  sinkholes  which  were  studied  were  selected  because  of  their 
location  at  the  low  end  of  a  large  uvala  or  karst  depression  (they  therefore 
receive  maximum  runoff  from  the  surface  watershed),  their  convenient  access, 
their  representation  of  both  collapse  and  subsidence  forms,  and  their  proxim¬ 
ity  to  each  other  (600  ft). 

Methodology 

10.  A  double  ring  of  electrodes  was  placed  circumferentially  around  the 
sinkholes,  using  copper-clad  steel  grounding  rod  cut  into  2-ft  lengths  and 
5/8  in.  in  diameter.  These  electrodes  were  driven  into  the  ground  to  a  depth 
of  18  in.  At  the  Harris  sinkhole,  the  inner  ring  consisted  of  eight  elec¬ 
trodes,  50  ft  in  diameter,  and  the  outer  ring  consisted  of  16  electrodes 

150  ft  in  diameter.  At  the  Hunter  sinkhole  the  inner  ring  of  electrodes  con¬ 
sisted  of  16  electrodes  in  a  30-ft  diameter,  and  the  outer  ring  consisted  of 
16  electrodes  150  ft  in  diameter.  The  measuring  electrodes  were  numbered 
starting  with  electrode  1  as  the  northern  electrode  in  the  inner  ring  and 
increasing  to  8  or  16,  proceeding  clockwise,  for  the  Harris  or  Hunter  sink¬ 
holes,  respectively.  In  the  outer  ring  the  most  northern  electrode  was  num¬ 
bered  9  or  17  and  increased  to  number  24  or  32  proceeding  clockwise  for  the 
Harris  or  Hunter  sinkholes,  respectively. 

11.  Two  reference  electrodes  were  used.  The  first  (reference  elec¬ 
trode  1)  was  located  21.5  ft  due  west  of  electrode  21  of  the  Harris  sinkhole. 
The  second  (reference  electrode  2)  was  located  in  the  center  of  the  Harris 
sinkhole  approximately  19  ft  below  surface  elevation.  The  reference  elec¬ 
trodes  were  the  same  size  and  material  as  the  other  electrodes,  and  they  were 
also  installed  in  the  same  manner. 

12.  SP  measurements  were  obtained  by  using  a  high  impedance,  digital 
voltmeter.  The  meter  was  connected  between  one  of  the  reference  electrodes 
and  a  reel  of  wire.  One  member  of  the  measuring  team  laid  out  wire  from  the 
reel  to  one  of  the  numbered  electrodes  located  around  the  sinkhole.  The  wire 
was  then  clipped  onto  the  electrode,  thus  completing  the  measurement  half  of 
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the  circuit  (the  earth  being  the  other  half),  and  the  Individual  relayed  by 
voice  the  number  of  the  specific  electrode  to  the  person  reading  the  voltmeter. 
The  person  reading  the  voltmeter  recorded  the  measurement,  disconnected  the 
reference  electrode,  and  took  a  second  measurement  using  the  second  reference 
electrode.  Upon  completion  of  this  procedure,  the  Instrument  reader  relayed 
by  voice  to  the  reel  person  who  then  proceeded  to  the  next  numbered  electrode 
where  the  measurement  procedure  was  repeated.  This  method  was  continued  until 
all  electrode  readings  using  both  reference  electrodes  were  taken. 

13.  Data  were  collected  approximately  biweekly  for  a  6  month  period. 

When  SP  data  around  a  sinkhole  revealed  a  consistent  and  repetitious  pattern, 
a  tracking  array  of  electrodes  was  Installed  beyond  the  edge  of  the  sinkhole 
to  determine  If  the  Inferred  flow  patterns  could  be  detected  at  some  distance 
away.  The  lowest  SP  readings  were  Inferred  to  be  subterranean  drainage  paths. 
The  tracking  array  consisted  of  a  grid  of  16  electrodes  spaced  on  50  ft  cen¬ 
ters.  The  16  electrodes  In  a  tracking  array  were  alphanumerlcally  labeled  A1 
through  PI  In  the  array's  first  position  and  A2  through  P2  In  the  array's  sec¬ 
ond  position. 

14.  Each  array  consisted  of  four  rows  of  four  electrodes.  The  north¬ 
ernmost  row  for  tracking  array  one  was  labeled  from  left  to  right  A1  through 
Dl.  The  other  rows  were  labeled  consistent  to  this  pattern  with  the  south¬ 
ernmost  row  having  electrodes  labeled  from  left  to  right  Ml  through  PI.  As 
the  array  was  advanced  to  a  second  position,  the  same  numbering  pattern  was 
followed.  Once  a  repetitious  pattern  was  discerned,  the  array  was  advanced  In 
the  preferred  direction  and  the  measurements  repeated.  In  this  way  It 
appeared  that  a  subsurface  flow  path  between  the  Harris  and  Hunter  sinks  was 
discernible . 

15.  In  addition  to  SP  data,  electrical  resistivity  measurements  were 
taken  at  the  site  for  the  purpose  of  providing  some  confirmation  of  the  SP 
results.  Using  the  already  Implanted  electrodes,  azimuthal  measurements  were 
obtained  by  rotating  a  50  ft  spacing  Wenner  electrode  array  22.5  deg  until 
eight  readings  were  taken.  This  rotational  measurement  scheme  resulted  In  a 
characterization  of  the  resistivity  variation  as  a  function  of  azimuth  around 
the  sinkhole  and  reflected  the  electrical  anlsotrophy  of  the  underlying  bed¬ 
rock  and  regolith. 

16.  Since  the  study  area  Is  located  In  a  very  rural  area,  no  studies 
were  conducted  for  stray  currents  or  interference.  Most  readings  were  very 


stable.  Occasional  erratic  readings  were  a  result  of  cattle  rubbing  on  the 
electrode  and  loosening  it  in  the  soil .  Further  tamping  corrected  these  prob¬ 
lems.  Air  temperature  and  precipitation  data  from  an  official  station  at 
Lexington  (6  miles  from  study  area)  are  used  to  evaluate  the  effects  of  these 
parameters  on  the  SP  readings. 

Results  and  Discussion 

17.  The  SP  data  from  both  sinkholes  that  were  studied  are  given  in 
Tables  1,  2,  and  3.  Table  4  lists  the  SP  values  for  electrodes  in  the  two 
tracking  arrays  that  were  used  between  the  sinkholes.  The  SP  data  for  a  given 
electrode  varied  greatly  over  time,  but  variation  between  electrodes  remained 
relatively  constant.  Changes  in  ground-water  flow  as  well  as  surface  ground 
temperature  and  precipitation  or  a  combination  of  these  and  other  factors  are 
probably  responsible  for  fluctuating  SP  values. 

18.  Since  the  raw  data  shown  in  Tables  1,  2,  3,  and  4  indicate  vari¬ 
ability  over  the  collection  period,  an  effort  was  made  to  filter  the  data  by 
using  ranks  to  preserve  the  ability  of  the  data  to  document  consistent  nega¬ 
tive  anomalies  and  to  remove  seemingly  erratic  swing  in  the  raw  data.  Nega¬ 
tive  SP  readings  for  a  given  collection  day  would  be  given  a  low  ranking,  e.g. 
the  electrode  having  the  most  negative  value  for  the  day  would  be  given  a  rank 
of  one,  and  the  electrode  having  the  most  positive  value  would  be  given  a  rank 
of  24  at  the  Harris  sinkhole  or  32  at  the  Hunter  sinkhole. 

19.  An  example  of  the  ranked  data  and  the  average  value  and  standard 
deviation  of  the  rank  data  over  the  collection  period  for  the  Harris  sinkhole 
using  reference  electrode  1  is  presented  in  Table  5.  Negative  readings  are 
shown  as  low  rankings  and  positive  readings  are  shown  as  high  rankings.  The 
standard  deviation  of  the  ranking  shows  the  relative  sensitivity  of  an  elec¬ 
trode  to  time  varying  effects  such  as  water  flow  while  systematic  changes  that 
are  believed  to  be  associated  with  temperature  effects  are  filtered  out  by  the 
ranking  procedure. 

20.  The  ranking  procedure  has  both  advantages  and  limitations.  The 
major  advantage  is  the  removal  of  systematic  changes  in  the  data.  These  sys¬ 
tematic  changes  translate  the  origin  for  the  data.  The  major  disadvantage  is 
the  Insensitivity  to  factors  which  exaggerate  scale.  Increased  ground-water 
flow  may  make  the  most  negative  readings  more  negative.  This  change  is  not 
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Table  3 

Values  of  Spontaneous  Potential  in  Millivolts  Using  Reference 


Electrodes  1  and  2  and  the  Circumferential  Positioned 


Electrodes  Around  the  Hunter  Sinkhole 


Oct  22 


Reference  1  Reference  2 


Dec  18 


Reference  1 


Jan  1 


Reference  1  Reference  2 
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visible  In  the  rankings.  Only  when  readings  (raw  data)  pass  each  other  do 
changes  In  ranking  occur. 

21.  Figure  2  shows  the  electrode  pattern  and  the  average  ranked  SP  val 
ues  at  the  Harris  sinkhole  using  reference  electrode  1.  A  preferred  drainage 
pattern,  based  on  selection  of  electrodes  with  the  lowest  ranked  values,  Is 
Indicated  In  Figure  2.  Figure  3  Is  an  enlarged  topographic  map  of  the  study 
area  which  Incorporates  data  from  both  sinkholes,  the  tracking  arrays  between 
them,  and  the  Inferred  subsurface  flow  path. 


SCALE 

25  0  25  50  75  FT 
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Figure  2.  Electrode  pattern  spacing,  average  ranked  spontaneous 
potential  values,  and  the  Inferred  flow  pattern  at  the  Harris 
sinkhole  using  reference  electrode  1 


Figure  3.  Topographic  map  of  the  Harris  and  Hunter  sinkholes 
and  their  associated  circumferential  ring  electrodes,  tracking 
arrays,  and  the  drainage  path  as  inferred  by  SP  data 


22.  The  inferred  subsurface  flow  path  connecting  the  Harris  and  Hunter 
sinkholes  Is  over  600  ft  In  length.  It  does  not  follow  the  direction  of  sur¬ 
face  drainage  but  appears  to  be  oriented  along  strike  of  the  underlying  bed¬ 
rock  (approximately  N  50  deg  E) .  This  strike  direction  also  appears  to 
control  passage  development  of  nearby  caves  as  well  as  sinkhole  distribution. 
Most  of  the  sinkholes  In  the  general  area  are  aligned  In  this  preferred  orien¬ 
tation  (see  Figure  1).  The  Hunter  sinkhole  Is  located  about  20  ft  higher  than 
the  Harris  sinkhole,  and  from  observation  during  runoff  events  It  apparently 
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Intercepts  surface  drainage  and  conducts  the  water  subterraneanly  (via  strike 
controlled  fissures)  into  the  lower  Harris  sinkhole.  SP  data  appear  to  con¬ 
firm  this  model. 

23.  Further  confirmation  of  the  SP  data  was  obtained  from  observation 
of  surface  runoff  that  is  channeled  by  a  road  and  cattle  guard  into  the  Harris 
sinkhole  during  periods  of  heavy  precipitation.  SP  data  (from  electrode  14 

to  electrode  2)  detected  this  flow  path  long  after  surface  runoff  had  ceased 

flowing.  This  flow  path  is  Inferred  to  be  a  recharge  path  within  the  regollth 

as  well  as  a  surface  runoff  path. 

24.  Table  6  lists  the  averages  and  standard  deviations  of  the  ranked  SP 

data  at  the  Harris  sinkhole  for  reference  electrodes  1  and  2.  From  Table  6  it 
is  evident  that  low  ranked  measuring  electrodes,  i.e.,  those  electrodes  having 
low  or  negative  SP  values,  are  generally  consistent  even  though  different  ref¬ 
erence  electrodes  are  being  used.  Actual  differences  in  SP  values  may  be  due 
to  the  elevation  difference  of  32  ft  between  reference  electrodes  (reference 
electrode  2  being  at  the  bottom  of  the  aven  in  soil  which  had  recently  col¬ 
lapsed  into  the  bottom  of  the  sinkhole) .  Differences  in  soil  moisture  or 
structure  of  the  in  situ  soil  at  reference  electrode  1  and  the  collapsed  soil 

at  reference  electrode  2  may  also  explain  the  range  of  SP  values. 
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25.  Azimuthal  electrical  resistivity  data  from  the  Harris  sinkhole  is 
shown  in  Figure  4.  Although  this  represents  only  one  set  of  data,  it  is 
Interesting  to  note  that  the  major  axis  of  the  elliptical  plot  reflects  the 
strike  direction.  Taylor  (1985)*  has  shown  that  joint  orientation  is  one  of 
the  most  useful  and  reliable  parameters  that  can  be  obtained  from  azimuthal 
resistivity  measurements.  The  character  of  this  ellipse  may  be  Interpreted  to 
determine  the  transverse  and  longitudinal  resistivity  for  the  subsurface  mate¬ 
rial  at  the  site.  Dominant  strike  controlled  jointing  and  secondary  dip  con¬ 
trol  (bedding  plane)  of  ground-water  movement  are  suggested  by  the  resistivity 
data  as  well  as  the  SP  data  and  observation  of  outcrops  and  cave  passages  in 
the  area. 

26.  Figure  5  relates  precipitation,  average  dally  air  temperature,  and 
SP  values  from  a  randomly  chosen  electrode.  The  qualitative  trends  of  SP  data 

*  R.  W.  Taylor.  1985.  "The  Determination  of  Joint  Orientation  and  Porosity 
from  Azimuthal  Resistivity  Measurement,"  a  paper  presented  at  the  Annual 
Meeting  of  the  Association  of  Engineering  Geologists,  Raleigh,  N.C. 
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Table  6 

The  Averages  and  Standard  Deviations  of  the  Ranked 
Spontaneous  Potential  Data  at  the  Harris  Sinkhole 
for  Reference  Electrodes  1  and  2 


Reference  Electrode  1 

Reference 

Electrode  2 

Average 

Standard 

Average 

Standard 

EL  No. 

Ranking 

Deviation 

Ranking 

Deviation 

1 

19.2 

3.6 

19.4 

0.4 

2 

10.4 

7.7 

12.1 

7.3 

3 

10.5 

2.8 

11.4 

2.9 

4 

12.1 

5.4 

10.7 

5.7 

5 

13.8 

5.8 

14.2 

5.8 

3.4 


3.2 


2.5 


IHI 


3000  0HM-FT 


2000  0HM-FT 


Figure  4.  Azimuthal  resistivity  results  in  ohm-ft 
taken  at  the  Harris  sinkhole  site 


for  other  electrodes  are  similar.  Temperature  and  precipitation  data  are  from 
an  official  weather  station  in  Lexington,  6  miles  from  the  study  area.  Pre¬ 
cipitation  values  are  total  5-day  accumulations  up  to  the  date  on  which  SP 
data  were  taken.  The  monthly  accumulations  are  represented  by  horizontal 
lines.  Average  daily  air  temperature  is  assumed  to  give  an  indication  of  sur¬ 
face  ground  temperature.  In  general,  it  seems  that  SP  values  decrease  with  a 
decrease  in  temperature  or  an  Increase  in  precipitation.  This  finding  is 
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Figure  5.  Precipitation,  air  temperature,  and  self  potential  value 
for  electrode  17  taken  over  the  measurement  dates  conducted  during 

this  study 

consistent  with  work  presented  by  other  investigators  in  that  an  Increase  in 
water  flow  resulted  in  a  lowering  of  SP  values.  The  effects  of  temperature 
fluctuation  on  SP  values  have  been  studied  by  Corwin  and  Hoover  (1978).  They 
have  noted  a  coupling  effect  and  describe  a  decrease  in  SP  values  over  geo¬ 
thermal  areas  (where  temperature  extremes  are  high).  In  this  study  area,  as 
surface  temperatures  decrease  in  winter  months,  a  thermal  couple  develops 
between  the  cool  surface  and  warmer  ground  at  depth.  This  thermal  difference 
may  approximate  the  effects  of  geothermal  areas,  and  thus  SP  values  decrease 
as  surface  temperature  decreases.  In  summer  months  the  temperature  differ¬ 
ences  between  surface  and  subsurface  are  reversed.  A  similar  correspondence 
between  temperature  and  SP  is  noted,  although  the  response  is  not  as  marked; 
this  could  probably  be  due  to  a  greater  lag  time  in  ground  warming  for  these 
months . 
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Conclusions 


The  major  conclusions  from  this  research  are  the  following: 

a.  The  SP  technique  used  in  this  study  was  effective  in  detecting 
SP  anomalies  which  were  inferred  to  represent  surface  and  sub¬ 
surface  drainage  or  flow  paths  around  a  sinkhole.  In  addition, 
the  SP  technique  was  able  to  track  one  anomaly  (flow  path)  for 
over  600  ft  into  another  sinkhole. 

b.  Confirmation  of  the  interpretation  of  the  SP  data  was  obtained 
by  visual  observation  of  surface  runoff  into  the  instrumented 
sinkhole,  by  electrical  resistivity  measurements,  and  by  geo¬ 
logical  studies. 

c.  Although  SP  values  varied  significantly  during  the  6-month 
(July  to  December  1985)  testing  period  for  individual  elec¬ 
trodes,  the  ranking  values  between  electrodes  were  consistent. 
This  was  also  true  if  the  position  of  the  reference  electrode 
was  changed.  It  appeared  that  changes  in  precipitation  and 
surface  ground  temperature  had  the  greatest  affect  on  the  vari¬ 
ation  of  SP  data  over  the  testing  period. 


Future  Studies 

28.  The  results  of  this  preliminary,  proof  of  concept  research  for 
studying  flow  patterns  around  sinkholes  and  for  geotechnical  applications  of 
the  SP  method  in  general  will  be  greatly  enhanced  by  studies  to: 

a.  Confirm  and  quantify  the  effects  of  major  environmental  factors 
on  SP  measurements  by  monitoring  electrode  arrays  over  long 
periods  of  time  (>  6  months). 

b.  Instrument  and  monitor  SP  electrodes  around  sinkholes  under  the 
controlled  release  of  known  quantities  of  water  and  at  a  sink¬ 
ing  stream  site  with  known  flow  rate  in  attempts  to  relate  SP 
response  to  flow  rate. 

c.  Compare  and  quantify  differences  in  the  response  of  metallic 
and  non-polarizing  electrodes  to  environmental  factors  and 
assess  the  utility  of  the  different  electrode  types  for  long 
term  monitoring  applications. 

d.  Develop  a  library  of  SP  modeling  tools  (computer  programs)  for 
predicting  SP  response  and  for  interpreting  field  survey  data. 

Items  a,  b,  and  c  will  be  addressed  in  Report  2  in  this  series,  and  items  a, 
c,  and  d  will  be  the  subject  of  Report  3  in  this  series. 
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Synopsis  of  Recommended  Field  Procedure 

29.  It  is  clear  from  the  proceeding  section  that  field  procedures  for 
geotechnical  applications  of  the  SP  method  are  still  evolving.  Environmental 
factors  such  as  vegetation,  cultural  features,  and  variations  in  temperature, 
soil  type,  and  moisture  content  over  the  survey  site  and  during  the  period  of 
a  long-term  monitoring  effort  can  apparently  effect  metal  electrodes,  such  as 
the  copper  clad  steel  rods  used  in  the  present  study,  to  a  much  greater  extent 
than  non-polarizing  electrodes.  While  it  may  be  possible  to  compensate  for 
the  environmental  effects  on  metal  electrodes  by  careful  field  procedures  and 
data  processing  techniques  such  as  presented  by  Butler  et  al.  (1984),  the 
extra  time  and  effort  involved  may  obviate  the  advantages  of  lower  cost,  easy 
availability  and  simplicity  of  installation  and  maintenance  of  the  metal  elec¬ 
trodes  compared  to  non-polarizing  electrodes.  With  due  consideration  to  the 
above  reservations  regarding  the  use  of  metal  electrodes,  the  following  is  an 
abbreviated  field  procedure  recommendation,  with  a  more  detailed  and  defini¬ 
tive  set  of  recommendations  to  be  given  in  Reports  2  and  3  of  this  series: 

a.  Plan  the  SP  survey  to  use  a  single  fixed  reference  electrode  if 
at  all  possible. 

b.  Use  electrodes  which  have  constant  length  and  install  (drive) 
each  to  a  constant  depth. 

c.  Install  the  electrode  array  at  least  48  hours  before  the  first 
set  of  SP  measurements  are  acquired. 

d.  Carefully  record  any  variations  in  soil  type  along  the  survey 
lines  and  the  locations  of  any  cultural  features  such  as 
piezometer  pipes  or  buried  pipes  or  conduits. 

e.  Note  any  unstable  readings  during  SP  surveys. 

f.  Keep  a  record  of  rainfall  times  and  amounts  during  the  course 
of  an  SP  monitoring  study. 

g.  Examine  each  set  of  SP  measurements  for  correlations  with  known 
conditions  and  for  evidence  of  reference  level  shifts  compared 
to  previous  survey  results. 

h.  Use  the  procedure  presented  in  Butler  et  al.  (1984)  to  process 
SP  monitoring  data  to  compensate  for  reference  level  shifts  and 
variations  in  soil  type,  moisture  content,  and  other  non¬ 
seepage  related  SP  variations;  this  procedure  produces  an  SP 
anomaly  data  set. 

1.  Correlate  the  SP  anomaly  data  with  reservoir  pool  level. 
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